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ABSTRACT 

Collective mode of one- dimensional Fermi gas has been 
evaluated over the whole coupling regime from weak coupling 
(high density) to strong coupling (low density). The evaluation has 
been performed using theoretical model of Alam and Schuck 
taking attractive 5 interaction potential under RPA approximation 
at zero temperature. We observe that in the weak-coupling limit 
particle-hole RPA approaches to low momenta and the collective 
mode in the strong coupling limit reproduces the Bogoliubov 
mode for the weakly interacting bosons. Our theoretically 
evaluated results are in good agreement with other theoretical 
workers. 

Keywords: collective mode, interaction potential, Random Phase 
Approximation (RPA) Quasi-particle, weak- coupling, strong 
coupling, Bogoliubov mode. 



INTRODUCTION 

The one dimensional Fermi gas with 
attractive 8 interaction among the fermions 
has been used as a model to know the 
physical properties of realistic Fermi system. 
The exact solution for its ground state 
energy is known from the Bathe ansatz 1 , 
therefore one can test approximate solutions 
of the systems. Quick, Esbagg and deLalno 2 



have developed mean field approximation 
basically of two types (a) plane wave Hartre- 
Fock (HF) (b) non-plane wave (HF) and 
BCS. It was found that the BCS solution can 
describe the cross-over between weak 
coupling (as a weakly interacting gas of 
fermions) and strong coupling (gas of 
bosonic two particle pairs) in the system. In 
particular, they find that the BCS solution 
for the ground state energy coincides with 
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the exact solution in both weak and strong 
coupling. In this sense, the system may serve 
as a simple model to study the transition 
between weak and strong coupling 
superconductivity in ID Fermi system. This 
transition between weak and strong coupling 
has been discovered by Leggett 3 and by 
Noziers and Schmitt Rink 4 in three 
dimensions. The simple form of the 
interaction allows one to carry out 
approximation beyond the mean-field level 
such as ordinary random phase 
approximation (RPA's) or generalized RPA's 
in controlled fashion. In this approximation, 
one is able to calculate contribution to the 
ground state energy of the system beyond 
the mean field ground state energy. Williams 
and Bloch 5 have also discussed the ordinary 
(particle hole) RPA for the one dimensional 
electron gas. The same approach was also 
applied by Bernner and Haug 6 . Friesen and 
Bergersen 7 have applied Singwi Sjolander 
generalization of RPA 8 to the ID electron 
gas. This model can serve as an initial 
approximation to realistic systems such as 
quasi one dimensional metal 9 . 

In this paper, using the theoretical 
formalism of T. Aim and P. Schuck 10 , we 
have studied the collective modes of the one 
dimensional Fermi gas within the quasi 
particle random phase approximation (RPA). 
One calculates the collective excitations for 
the ID attractive Fermi gas with a 5 
interaction at T = 0 by applying the quasi 
particle RPA. One obtains quasi particle 
RPA equation for the two quasi particle 
propagator starting from the Bogoliubov 
transformed Hamiltonian of the system 
using the equation of motion for the two 
particle Green's function. The homogeneous 
two particle equation yields the condition for 



the collective excitation in the system. It 
coincides with result found by Anderson 11 , 
Rickayzen 12 , Bardasis and Schrieffer 13 and 
others 14 ' 15 by the equation of motion method. 

MATHEMATICAL FORMULAE USED 
IN THE STUDY 

One writes the Hamiltonian for one 
dimensional Fermi gas in second 
quantization 

h = I> (k)c + k c- k +1 £<^#i^oqc; 2 Q C t< 

(1) 

where k denote momentum and spin 
quantum number of the particle and 
{k l k 2 \V \K } K 4 ) j s ^ anti-symmetrized 
matrix element of the two body interaction. 
Using Bogoliubov transformation of the 
creation and annihilation operators 
Hamiltonian (1) is transformed into a new 
form. The transformed Hamiltonian is 
derived by several author, 16-18 The 
transformed Hamiltonian is written in the 
following form 

H = H" + Y J HlX < a ki + i £ (H™ 2 a\ a + k2 + H. c) 

k,k 2 1 k,k 2 

k,k 2 k 2 k t 

+ ^ Tj H hk 2 k,k, ^k^kl^k^k, 

(2) 

In equation (2), H° is the BCS 
ground state energy and H 11 is the diagonal 
part, H 20 is the off-diagonal part of the 
Hamiltonian. If one demand that the off- 
diagonal part H 20 vanishes then one obtains 
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the well known relations of BCS gap 

equation 19 . G l Q l (kq,a))= "' 

The H 31 term in the Hamiltonian in 
eq n (2) does not contribute to the RPA 
equations. The other terms H 40 and H 22 in 
the Hamiltonian describe the residual 
interaction among the quasi particle. 

These terms are neglected in the w h ere 
model BCS approximation. Therefore in 
order to go beyond the BCS mean field 
approximation one has to include the 
residual interaction among the quasi particle. 
Now, one can treat the quasi particle within expressions 
the generalized RPA approximation. For G u (kk'q,co) + G 21 (kk'q,co) 
this, one introduces two particles Green's j 
function with respect to quasi particle basis = G 1 0 1 (kq,G))S kk ,+j— 2 —^2E kq m(k,q)\^ rq 



Ek, q ~ E k + E k+q 

We have assumed that collective 
pairs have zero total spin One arrives on the 



and derive the equation of motion for the 
Green's function G. The equation has the 
form of Dyson equation 20 for the two 
particle propagator matrix G. The static anc j 
elements of the mass operator A and B are 
given by the double commutation 21 



+E^n{k,q)K +ad{k,qf VqA 



(6) 



G n (kk'q,co)-G 2l (kk'q,(o) 
1 



Gl\kq,Q))S kX +- 



[2(m{k,q)]Z k , qia 



k\k 2 k 3 k 4 

and 

^k } k 2 k 3 k 4 



(3) 



(4) 



Here > ' refers to averaging with respect to 
the BCS ground state. But in place of BCS 
ground state one has correlated ground state 
which corresponds to a generalized quasi 
particle RPA 



22 



+G*l{k,q)A k , w +E k j{k,q)r k ^ a) 

(7) 

The quantities Z, A an d T are given by the 
equations 

Z r q , a =~I H k "4° n ( k " k <l> W ) + G2> ( k " k «> a )] 

z r . 

(8a) 

A^^-v^nik^qlG^k-kq^+G^k-kq^)] 

k" 

(8b) 

r*„ = ~ v % l)[G l \k"kq, co) + G 2 ' (k"kq, co)] 

(8c) 



Taking 5 interaction one solves the The quantities m(k,q), n(k,q) and l(k,q) are 
system equation in the (0 representation. One combination of u k and v k and is given by 



obtains coupling equation for G 11 and G 21 
also for G 22 and G 21 . One obtains 



m(k,q) = u k v k+q +v k u k+q 
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(9b) 



n(k,q) = u k u k+q +v k v k+q 



£(k,q)-- 
and 



U k U k+q V k V k+q 



1-v 



E„ 



Where 



E k =je 2 k +A 2 k 



(9c) 



(9d) 



(9e) 



Multiplying in equation (7) 
subsequently with m with n and in equation 

(7) with ^ and summing over k, one arrives 
at the system of equation for the quantities 
Z, A and r 



l+v/ E ,„,„(<7,ft>) 



v 2 
f 



\+vI Ett {q,G>) 2vl afm (q,0)) 
V -h,n,M><°) X + vI Emm (q,G>) 

A \ \-vn(k,q)Gi; ^ 
x T = -vl(k,q)Gl x 

-v/2m(k,q)G l Q l j 

The quantity I a b c are in the notation 



(10) 



_y a(k,q)b(k,q)c(k,q) 

a,b,c ~ Z-i I ,J1 t?2 \ 



(11) 



Here a(k,q) = [E kq (0] and b(k,q), c(k,q) = 

[n(k,q), ^ (k,q), m(k,q)]. This is a linear in 
homogenous system of equation for the 
quantities A, T and Z. This can be solved 



with the help of matrix inversion. Equation 
(10) is a eigen value problem for the 
determination of the collective modes in the 
quasi particle RPA. The condition for the 
nontrivial solution is the vanishing of 
determinant. 



+ vI EM (q,Sl) vl nilt (q,£l) vI Enm (q,Ci) 
vI nnl {q,Q.) l+vI w (q,Cl) 2vI ai Jq,0.) =0 
v/2/ £ ,„,„, (<?,£!) v/2I nnl {q,a) l + vI Emm (q,£l) 

(12) 

where Q. (q) denotes the eigen value for the 
collective excitation. 



Weak coupling case 

The analysis of the weak coupling 
case has been given by Belkker and 
Randeria 23 . One quotes their result for the 
weak coupling collective made in one 
dimension. 



Q(q) = 



1 m 
1-v 

71 k r 



1/2 



k l' 

—q 

m 



(13) 



where the density of states in one dimension 
for parabolic dispersion was used. The long 
wavelength collective modes in the weak 
coupling case have a phonon like spectrum 
and are independent of the gap. If one 
compares equation (13) with the small-q 
expansion of the collective modes in the 
particle hole RPA, one finds that both 
coincides. This means that in weak coupling 
the behavior of collective modes for small q 
is not changed from the normal particle hole 
RPA. 

However, for large q and in 
particular near the point q=2k F the quasi- 
particle differs from the particle hole RPA. 
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Strong coupling limit 

Now, from equation (12) one can 
expand the determinant of the collective 
modes for small -q and - Q. limit with 
respect to small gap. In the strong coupling 
limit the gap goes to zero with density. One 
expands the gap equation in terms of A 2 



1 = — \dk 

27T J 



1 



If. 



1 + 



1/2 



I6| , 



= — \dk-^— -v — f dk - = — /, — ^—J, 



(14) 



Now in equation (14) the integral in 
the first line is convergent for a contact 
interaction without a cutoff due to the one 
dimensionality of the system. An analogous 
expansion of the BCS density equation 
yields 



1 9 

n = — AX, 
2n 



(16) 



Now, the effective chemical potential can be 
expressed in terms of the density and the 
coupling strength as 



(^*) 1/2 =^ m v2 v± JA mv 2_^ v 



8 



3 

— i 
8 



(17) 



Here one has introduced integral J ; . It has 
been evaluated for the strong coupling limit 
i.e. for \L *<0. Here |i * is effective chemical 
potential including the quasi particle shift 

°° 1 
J i =jdk=- 

0 



2m 



+ jU 



7 i 2 = J dk — 



2m 



+ JU 



Then 

J, = W /2 /(2|//|) 
J 2 =m 1/2 ^/(2|//*|) 
J 3 = 3m 1/2 x/(2 7/ Y\) 



1/2 
.|\ 3 / 2 



5/2 



(15a) 



(15b) 



(15c) 



In the limit of zero density or zero gap, 
respectively equation (17) yields the 
condition 



2\ju*\ = 



mv 



= E n 



(18) 



This shows that in the extreme strong 
coupling limit the chemical potential i.e. the 
energy to remove a particle from the system 
in just half the two particle binding energy - 
E 0 in the Vacuum. 2 

The long wavelength dispersion relation in 
the strong coupling limit is given by 

<%) = cq (19a) 

where the sound velocity c depends on the 
gap A. Now using the low density expansion 
of the BCS density (16) to substitute A 2 by 
the density n, one obtain the expression for 
the collective modes in the strong coupling 
limit. 
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vn 
4m 



1/2 



a{ q ). 

(19b) 

Introducing the pair mass m B = 2m and the 
pair density n B = n/2 the above expression 
reduces to 

r n 1 / 2 
\ vn R 

a{q)= ^ q 

L J (19c) 

Equation 19(c) is the well known 

Bogoliubov dispersion relation for the 

weakly interacting Bose gas ' 5 in the limit 

of small q which is linear in q i.e. phonon 

like. Thus starting from interacting fermions 

with an attractive interaction the quasi 

particle RPA in the strong coupling limit 

yields the dispersion relation for weakly 

interesting gas of bosons (two particle bound 

states). The magnitude of the repulsive 

interaction among the bosons in equation 

19(c) is given by the fermionic interaction 

strength v and is consistent with the result of 

Hanssmann 26 and others. 27 32 

Table T t 

An Evauated results of collective mode Q. for 
small momenta q in the weak coupling limit 
(u*=3.386 E 0 ) The numerical solution is 
compared with weak coupling expansion (13). 



q/k F 


QJ 


|i 


Numerical solution 


Weak coupling 




Equation (12) 


expansion eq. (13) 


0 


0.00 


0.00 


0.1 


0.126 


0.132 


0.2 


0.268 


0.275 


0.3 


0.355 


0.368 


0.4 


0.568 


0.599 


0.5 


0.573 


0.702 


0.6 


0.626 


0.925 


0.65 


0.645 


1.027 


0.70 


0.689 


1.146 


0.75 


0.735 


1.207 


0.80 


0.794 


1.258 



DISCUSSION OF RESULTS 

In this paper, using the theoretical 
formalism of T. Aim and P. Schuck, 10 we 
have studied the behavior of collective 
modes over the coupling regime weak 
coupling (high density) to strong coupling 
(low density). We observe that the treatment 
of the residual interaction in the Hamiltonian 
(2) within the quasi particle or generalized 
RPA allows one to study the behavior of the 
collective over the whole coupling range. 
Collective modes in one dimension were 
found for the case of an attractive 8 
interaction. One observes that in the weak 
coupling limit one recovers Anderson's 
results 11 whereas in the strong coupling limit 
the Bogoliubov dispersion relation 24 for the 
interacting Bose gas of two particle pairs can 
be derived from the quasi particle RPA. This 
is consistent with the fact that the BCS 
theory is capable of describing the extreme 
strong coupling limit i.e. the gas of two 
particle bound state properly and reproduces 
the exact result for the ground state energy 
in this limit. 2 In table T{ we have presented 
for small momenta q in the weak coupling 
limit (u*=3.386E 0 ). The numerical solution 
equation (12) is compared with weak 
coupling expansion equation (13). Both 
solution coincides at q/kF = 0.4. This shows 
the consistency of the numerical solution 
with the well known weak coupling result 
which was obtained by Anderson 11 in the 3D 
case. In table T 2 , we have presented the 
evaluated results of collective mode Q. for 
small momentum q in the strong coupling 
case (u* = (-) 0.4966 E 0 )' The numerical 
solution equation (12) is compared with the 
strong coupling expansion (equation 19c) 
and also with free particle dispersion Q. 
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=q 2 /2m B . We observe that the numerical 
solution starts linearly in q and is consistent 
with the strong coupling expansion gives in 
equation 19(c). This confirms the 
interpretation of the collective excitations in 
the strong coupling limit as Bogoliubov 
sound modes of the two particle Bose gas 
that is formed in the limit. The free particle 
solution reached the full solution for large q. 
Some recent 33 37 calculations also reveal the 
same type of behavior. 

Table T 2 

An Evauated results of collective mode £2 for 
small momenta q in the strong coupling limit 
(u*=(-)0.4966 E 0 ) The numerical solution 
[eq.(12)] is compared to the strong coupling 
expansion [eq.(19c)] and free particle dispersion 
£2 =q 2 /2m B (a 0 = Bohr radius) 







q/a 0 


Numerical 
Solution 
Eq.(12) 


Strong 
expansion 
eq.(19c) 


Free particle 
dispersion 
£2 = q 2 /2m B 


0 


0.000 


0.000 


0.000 


0.02 


0.0022 


0.0020 


0.0002 


0.04 


0.0032 


0.0032 


0.0004 


0.05 


0.0034 


0.0035 


0.0005 


0.06 


0.0037 


0.00355 


0.00282 


0.07 


0.0038 


0.00362 


0.00294 


0.08 


0.0042 


0.00375 


0.00307 


0.09 


0.0048 


0.00382 


0.00322 


0.10 


0.0052 


0.00397 


0.00375 


0.11 


0.0056 


0.00416 


0.00400 


0.12 


0.0062 


0.00438 


0.00425 


0.15 


0.0094 


0.00469 


0.00473 


0.20 


0.0126 


0.00525 


0.00098 
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